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We propose using nonlinear Mach-Zehnder interferometer (NMZI) to efficiently prepare photonic 
quantum states from a classical input. We first analytically investigate the simple NMZI that 
can filtrate single photon state from weak coherent state by preferrentially blocking two-photon 
component. As a generalization, we show that the cascaded NMZI can deterministically extract 
arbitrary quantum state from a strong coherent state. Finally, we numerically demonstrate that the 
cascaded NMZI can be very efficient in both the input power and the level of cascade. The protocol 
of quantum state preparation with NMZI can be extended to various systems of bosonic modes. 


Introduction.- Integrated photonics can achieve un¬ 
precedented interferometric stability (l, 2] and build large 
scale interferometers Hi- However, reliable quantum 
state preparation for integrated photonics remains an im¬ 
portant challenge, because interferometers and coherent 
input states are insufficient for quantum states prepa¬ 
ration. We may use either post-selection or nonlinear 
interaction to overcome this challenge. The approach of 
post-selection only requires linear optical elements and 
photon detectors, but the preparation of quantum state 
is probabilistic and conditioned on the outcome of the 
projective measurement m 3 . The approach of nonlin¬ 
ear interaction assisted by an ancillary two-level system 
(TLS) can deterministically prepare arbitrary quantum 
state of the photonic mode [8H10j. but it requires strong 
coupling between the optical mode with the single TLS, 
which is experimentally challenging for integrated pho¬ 
tonics. Alternatively, we may consider using the nonlin¬ 
ear optical waveguide combined with ultra-stable inter¬ 
ferometers to achieve reliable quantum state preparation, 
without requiring TLS post-selection [TlHl3j |. nor 

feedback/feedforward control [141] . 


In this Letter, we propose to use interferometry com¬ 
bined with Kerr nonlinearity to filtrate single photons or 
extract any desired quantum states from coherent state 
input, as illustrated in Fig. [2(a). We first present the 
idea of quantum state filtration (QSF) of single pho¬ 
tons, which keeps the desired single photon component 
by blocking the undesired component to a different port. 
We then generalize the idea to quantum state extraction 
(QSE), which not only keeps the desired component, but 
also extracts the desired component from the undesired 
component before blocking/redirecting the residual pho¬ 
tons. 


Single photon filtration.- We first consider the simple 
task of QSE of single photons. As shown in Fig. QJb), we 
use a nonlinear Mach-Zehnder interferometer (NMZI), 
with a Kerr nonlinear medium in one of the arms. Since 
Kerr nonlinearity can induce photon number dependent 
phase shift, we can design the NMZI to induce destruc- 
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Figure 1: (color online) (a) Schematic illustration the arbi¬ 
trary quantum state filtration and extraction from coherent 
state input, (b) The configuration for QSF of single photon 
from coherent state input |a), using the simple NMZI (con¬ 
sisting of Mach-Zehnder interferometer, and Kerr medium 
and phase shifter), (c) Three processes for two photon out¬ 
put of Path A for weak coherent input, (d) The probabilities 
of n photons output of Path A against the phase difference 
between two arms 4 >, with = 0.1 and a = 0.1. (e) The 
second-order correlation function ( g (2 ^) of light output of Path 
A against <f> for various tp with a = 0.1. 


five interference at the output port when there are two 
photons. More specifically, with a vacuum input at Path 
A (upper path) and a coherent state input at Path B 
(lower path), the input state to the filtration is 

Win = \vac) A ® \a) B , (1) 

where |a) = e _ l“l 2 / 2 ^r(^) n | vac), and a(a t ) and 
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fr(frt) are annihilation (creation) operators for Paths A 
and B, respectively. Each beam splitter (BS) induces a 
unitary evolution, 

Ubs{ 9 12 )= e Wl - 2iatb+abt \ (2) 

with 6 \ and d 2 for BS 1 and BS 2 , respectively. The evo¬ 
lution in the nonlinear Kerr medium in Path A is 

U K (fi, p) = e^“ ta+i¥3atataa , ( 3 ) 


where p is the Kerr coefficient and <f> is linear phase shift 
(relative to Path B). The final output state of the single 
photon filtration is 


Wont =Ubs{9 2 )U K ((I), <p)U B s(9i) W in 

OO OO 

= Y^^ 2 ^p,qW) P( W) q \ va c) , ( 4 ) 

p—0 q=0 

With = ELo Y!n~l Kp+q-n (?) ( P p- 7 ") X 

{-l) n - l (sm 0 2 ) p+n - 2 l (cose 2 ) q - n+21 . The probability 
of p photons at the output of Path A is 


P p = (p\ Tr s {| ip 2 ) W\} | p) = Wa? ■ 

q—0 

The second-order correlation function (T^] is 

(2 ) (a f a f aa) E^=2 p(p ~ x ) x P p ,^ 

g • (5) 

which characterizes the generated single photon state. 

For a weak coherent input \a \ 2 <C 1 , we have P 2 -C Pi 
and can safely neglect the probability of multiple photons 
(Pn> 3 )- By considering the leading contribution, we have 
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( 2 ) 


2P 2 „ 

2/72,0 

2 

1 - e i2 ^ 

Pi 2 ~ 

M 1,0 


(1 — rye - *^) 2 


( 6 ) 


where /ii,o = acosdicos02(e^ — rf) and p 2,0 = 
i(acos 0 icos 0 2 ) 2 [— + rj 2 ] for the simple 
NMZI with rj = tanditandg- The three terms in p,2,o cor¬ 
respond to three different processes with two photons at 
the output of Path A, as shown in Fig. Etc). The inter¬ 
ference of these three processes can be controlled by the 
linear phase shift (cf>) and nonlinear coefficient (p). The 
optimal condition for p 2 Q = 0 is 

= 1 ± \/l — e i 2 , P, (7) 

which can always be fulfilled as long as p 7^ 0, so that 
the leading contribution to g^ can be eliminated. 

Fig. 0 d) shows the probability of n photons at the out¬ 
put of Path A (P„) depending on the linear phase shift <j>, 
with parameters p = 0 . 1 , r\ = | 1 —Vl — e i2 ^| and a = 0 . 1 . 
We find that the P 2 is greatly suppressed for cf> ss 0 . 137 T, 


while the dominant single photon emission P\ P2,3,4 
is not significantly affected. In Fig. 0 e), the relation be¬ 
tween g( 2 ^ and tf> are plotted for different values of non¬ 
linear coefficient p, with a = 0.1 and 77 given by optimal 
condition from Eq. 0 . We find good agreement between 
the approximated analytical solution from Eqs. ©&0 
(solid lines) and the exact numerical solution from Eq. 
© (dashed lines). With increasing nonlinear coefficient 
ip, the deviation from g^ = 1 becomes more signifi¬ 
cant, due to the Fano interference of the three processes 
(Fig. 0 c)) contributing to p 2 ^. These Fano-like curves 
show sub-Poisson statistic with g^ « 0 for <j> close to 
the optimal condition (Eq. 0), where the two photon 
output can be totally forbidden due to destructive in¬ 
terference. Meanwhile, we can also find the constructive 
interference of the two-photon output, which gives rise to 
super-Poisson statistic (g^( 0 ) » 1 ) output. Comparing 
the curves with different nonlinear effect coefficients, the 
single photon filtration is more sensitive to phase <f> for 
smaller p, indicating the crucial role of nonlinearity. 

For QSF of single photon, the fidelity is F = Pi = 
(cosdicos02) 2 |ck 2 | |l — e l2lp \. The optimal condition re¬ 
quires 77 = tan0itan0 2 ~ 1, we have |cos0icos0 2 | < \ and 
Pi < p\a 2 \/ 2 , which implies that the fidelity depends on 
both the the Kerr nonlinearity coefficient and the inten¬ 
sity of the coherent state input. QSF with simple NMZI 
cannot suppress the components with n > 2 photons (see 
Fig. 0 d) and also Ref. 3 Z|), and it only works for weak 
coherent state |a 2 | -C 1, which significantly limits the 
fidelity. Moreover, the fidelity of QSF is fundamentally 
limited by the overlap between the input state and the 
target state, P succ < \ (ip 

out l^m)| 2 , because it blocks all 
undesired components. To go beyond this limit, we need 
to generalize QSF to QSE, which not only keeps the de¬ 
sired component, but also extracts the desired component 
from the undesired ones. 

Cascaded filtration.- To implement QSE, we consider 
the cascaded NMZI, with a series of NMZIs connected 
sequentially. As shown in Fig. 0 a), the basic element 
consists of a BS ( 0 ) followed by a linear phase shifter 
(</>) and a Kerr medium (p) in the upper path. The 
basic element can be represented by a standard two-port 
unitary [Fig. 0 b)] 


U{fi,p,0) = £7k(0, p)U B s- 


( 8 ) 


The cascaded MNZI with N elements can be character¬ 
ized by 

u N = n^L 1 u(cj> l ,pi,0i). ( 9 ) 

For example, the simple NMZI (Fig 0 b)) consists of N = 
2 basic elements, with (f > 2 = p 2 =0. 

The cascaded NMZI can not only keep the desired 
single-photon component, but also extract the (desired) 
single-photon state from (undesired) nmlti-photon states, 
as long as there are enough photons in the undesired com¬ 
ponent. We numerically optimize the fidelity by tuning 
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Figure 3: QSE for | faarget) = (|0) + \l))/y/2. (a) The prob¬ 
ability of |0) and 11), (b) the Fidelity ( F ) and Purity ( Q ) of 
the output. The results are obtained by optimize the param¬ 
eters of a chain of IV = 20 units under the constrain that 
1 - Po - Pi < 0.01. 


Figure 2: (color online) Cascaded NMZI. (a) The basic ele¬ 
ment consists of a BS ( 8 ) followed by a linear phase shifter 
(< t >) and a Kerr medium (i p) in the upper path, (b) Schematic 
representation of the element for cascaded NMZI. (c) Fidelity 
of the single photon extraction increases with the number of 
cascade elements, N. The parameters are optimized numeri¬ 
cally under the constrain P „>2 < 0.01 with ip = 0.1. (d) Fi¬ 
delity of Fock state extraction (n = 1, 2, 3) increase with |a| 2 
for cascaded NMZI with N = 40. The results are obtained 
by optimize the parameters of each unit under the constrain 
that 1 — Po — P„ < 0.01. 


the parameters of the N elements. As illustrated in Fig. 
He), the optimized fidelity of single photon extraction 
F = P\ increases with N monotonically, with asymptotic 
value F — > 1 — |(0|a)| 2 (dashed lines), because our pas¬ 
sive device cannot extract single photon from the vacuum 
component. Furthermore, the cascaded NMZI can ex¬ 
tract Fock state | n) with n = 1, 2,3, • • •. The asymptotic 
fidelity of n-photon extraction is F — > 1 — l( m l a )l > 

which can be achieved for \a\ 2 < 1.5 jn with cascaded 
NMZI of N = 40 elements, as shown in Fig. H^)- 

Arbitrary state Extraction. -Remarkably, the cascaded 
NMZI can extract arbitrary superposition of Fock 
states with a large coherent state input (|a| 1) with 

almost perfect fidelity. For 9i -C 1 with l = 1, ■ • • ,N, 
almost all input photons will be guided in Path B, 
which effectively remains as a coherent state (with small 
deviation of O ( 9)) for all intermediate stages. The effect 
of each beam splitter to the upper path can be regarded 
as an effective displacement operation to Path A, as 
D(ei) = e tia,] ~ tia with e/ = a9i and a small deviation 
of 0(e 2 /a 2 ) [l6|. In addition, the linear phase shift 
and Kerr nonlinearity can achieve the unitary evolution 
Uk{<Pi,v) = e l(/>iaia+i<pa,fafaa . Hence, the cascaded 
NMZI of N elements can induce the unitary evolution 


U K ( ct>N , ip) u (ejv) • • • Uk (fa, ip) u (e 2 ) U K (fa U (ei), 

which in principle can accomplish any desired unitary 
transformation for sufficiently large N and carefully cho¬ 
sen {fa, ez} i=1 ... N Uni- Despite the large overhead 


in N, this provides a generic approach using cascaded 
NMZI to extract arbitrary superposition of Fock states 
from a large coherent state with almost perfect fidelity. 

In practice, it is favorable to design the cascaded 
NMZI with a small number of elements. To illustrate 
the feasibility, we consider the target state | faarget) = 
(|0) + 11))/a/ 2 using N = 20 cascaded elements optimize 
the fidelity by tuning parameters of {tpi, fa, 9i } l=1 ... N . 
As illustrated in Fig. [3l we can improve the fidelity F 
and purity Q = Tr (p\) of the extracted state by increas¬ 
ing \a\ 2 . Both F and Q are greater than 97.5% when 
\a\ 2 > 1.5. It’s intriguing that a high fidelity QSE of 
superposition of Fock states can be achieved using a rea¬ 
sonable size coherent state and a finite-stage cascaded 
NMZI. 

Discussion.- The photonic integrated circuits (PICs) 
provides a promising platform for realizing cascaded 
NMZI, where arrays of beam splitters and phase shifters 
can be integrated on a chip [1, 0|. QSE can provide 
arbitrary input photon states for quantum information 
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processing yj, [7||. For experiment realization, the most 
challenging part is the Kerr nonlinear at single photon 
level (ip = 0.1 in this paper). Compared with schemes 
using single TLS at near-resonance condition, the Kerr 
nonlinearity can be realized by collective effect of ensem¬ 
bles and is more feasible for experiments. Recently, the 
interface between atomic ensemble and photonic waveg¬ 
uide for single photon nonlinearity, such as nanofiber 
0EI, hollow-core photonic crystal fiber [H, Q and 
integrated waveguide 24|, have been deonstrated [25- 
[27j . Novel mechanisms including giant Kerr effect [28| 
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and strong interaction between Rydberg atoms 1291 ] have 
also been proposed and demonstrated. As the integrated 
photonic chips starting the new era of fj-level (1000 pho¬ 
tons) nonlinear effect 0 53 , new materials, such as 
graphene 321 and topological insulator 53, might en¬ 
able next generation of PIC with nonlinearity at aJ-level 
( 10 photons). Therefore, the QSE is promising approach 
for future quantum state generation in PIC. 

The idea of QSE can be extended from optical fre¬ 
quency to microwave and terahertz frequencies. In par¬ 
ticular, the superconducting quantum circuits is read¬ 
ily to realize the QSE 0, 351. The mechanism of the 
QSE is very general, can also be generalized to other 
Bosonic excitations, such as surface plasmon, exciton- 
polariton, magnon and phonon. For example, quantum 
single phonon sources are possible by nonlinear mechan¬ 
ical oscillators (36]- 

Conclusion.- We have demonstrated that the simple 
NMZI can filtrate single photon state from a weak coher¬ 
ent state. Using cascaded NMZI, we can reliably extract 
arbitrary quantum state from a strong coherent state. 
Since our scheme only requires Kerr nonlinearity, linear 
phase shifter and beam splitter, it can be implemented 
in superconducting circuits, coupled optomechanical sys¬ 
tems, as well as photonic integrated circuits. 
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I. LIMITATION OF TWO-UNIT QUANTUM STATE FILTRATION 



Figure SI: The probability of n-photon (a) and the photon statistic g ^ (b) at output for increasing input coherent laser 
intensity \a 2 \. The solid lines are results for optimal conditions from Eq. 0, the dashed lines are numerically optimized results 
for constrains that P n >2 < 0.01 and P „>2 < 0.001. The dependence of Pi (c) and g ^ (d) against the optimization constrain 
P > 2. 


Since the fidelity of single photon filtration (the brightness of the single photon source) is limited as Pi < <p\ot 2 \/2, 
here we study the performance of the QSF with increasing input power \a 2 \ numerically. Shown in Fig. ISlI al by solid 
lines, the probabilities of Fock state outputs P n (n = 0,1, 2,3,4) against \a 2 \ is calculated, for Lp = 0.1 and other 
parameters satisfy the optimal condition Eq. 0. For weak input, the P\ increases linearly with |o:| 2 as expected, and 
g( 2 ) also increases [Fig. ISll bl], This means that the practical performance of the SPF deviated from the expected 
perfect single photon source that g^ ss 0 for \a 2 \ > 1. The imperfection for larger input coherent laser intensity 
should be attributed to the contributions of multiple photons with n > 3. This can be inferred from the probability 
of P 2 and P 3 in Fig.EHa), where the P 2 is approximately linearly depends on P 3 . This indicating that the 2-photon 
output comes from the 3-photon component of input state, which can’t be eliminated efficiently by two-unit QSF. 

Since the derivation of the optimal conditions is based on the assumption that \a 2 \ -C 1, the optimal condition 
may not be valid when increase the input coherent laser power. Therefore, to gain Pi as large as possible but keep 
the multiple photon probability as small as possible, we may optimize the parameters of QSF numerically for larger 
\a 2 \. For example, we optimize the Pi under the constrain P ra >2 < 0.01(0.001), the results are shown by dashed 
(doted) line in Figs. ISTT af and (b). For \a 2 \ > 1, the performance of the QSF in both Pi and g^ are improved. For 
example, with the input mean photon number \a\ 2 = 1, Pi = 0.047 and g^ = 0.48 for the optimal condition from Eq. 
0- With numerical optimization, we obtained the improved performance as Pl = 0.18 and g^ = 0.48 for constrain 
P n >2 < 0.01 and Pi = 0.10 and g^ = 0.22 for constrain P n >2 < 0.001. It’s obvious that, As illustrated in Fig. 2(c) 
and (d) by comparing the optimized Pl and </ 2 ' for different constrain P n > 2 , there is a trade-off relation between Pl 
and g 1 ' 2 ' 1 due to the lack of ability to suppress 2-photon output from multiphoton input state. 
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II. THE DEPENDENCE ON KERR NONLINEARITY 




(p N 

Figure S2: (a) The probability of single photon filtration by QSF against Kerr nonlinearity p for different input coherent laser 
intensity \a 2 \. (b) The probability of single photon extraction by TV-unit QSE against unit number N for different p. All results 
are optimized under constrains that P n >2 < 0.01. 


In the main body of the paper, most results of QSF are studied for fixed nonlinearity parameter p = 0.1. From the 
results of QSF, the fidelity of single photon filtration is also limited by the Kerr nonlinearity in the NMZI. Therefore, 
we provide more results to study the dependence of the performance of QSF/QSE on p. 

Shown in Fig. IS21 al are the fidelity of single photon filtration bu QSF against p for different \a 2 \. For both 
\a 2 \ = 0.5 and \a 2 \ = 1.0, the Pi increase with p, confirms that the asymptotic formula Pi < p\a 2 \/2, indicating 
that the performance of QSF by the simple NMZI can be improved by larger nonlinearity. By further increase p , Pi 
reaches the saturation value when p ss 1. 

In Fig. IS21 bl. the Pi by TV-unit QSE against N for various p are studied. It’s not surprising that the performance 
of QSE is also improved by increasing N, and reaches the saturation value 1 — |(0|cr}| . Comparing different curves, 
the higher the nonlinearity is, the less units required to reach the maximum fidelity. 



Figure S3: The P 2 for different input coherent laser power. 


III. IMPERFECTIONS 


To test the performances of QSE against the imperfections of parameters, we studied the performance of the single 
photon extraction with the optimized parameters for the QSE consist of N = 20 units with \a\ 2 = 1.0 and p = 0.1. 
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Figure S4: (Left) The Pi for different samples of varied parameters. (Right) The histogram of the Pi for the random variation 
of parameters. 


In Fig. S3, the P\ for varying input coherent laser intensity is shown. In a wide range of input coherent state mean 
photon number, the P\ is not deviated much from the optimized value. 

In Fig. S4, the Pi is tested with random perturbations of <f> and 6 from the optimized value. The perturbations 
are randomly and uniformly distributed in the range from — 0.01-7T to 0.017T. From the results, the performance of the 
QSE show certain tolerance to the parameter imperfections. 


IV. PROOF OF HIGH FIDELITY 


As evidences shown in Figs. (2)&(3), the purity and fidelity of the QSF/QSE can approach unitary when \a\ 2 and 
N approaches infinity. Here, we provide the analytical analysis to support this argument. In the example, the target 
state is 

| ^target) = (10) + |1))/V^. (S.l) 

In general, the output state can be represent by 

oo oo 

U |0) A |a ) B = |0) A c n (o,n \n) B + |1)^ c nCi, n _i \n — 1) B , (S-2) 

n=0 n =1 

where U is unitary transformation, c n = ^=e _ l“l 2 / 2 for coherent input jcu), and the coefficients satisfy|Co,„| + |Ci,n| = 1- 
For N 1 , it’s possible to achieve near-perfect extraction of Fock states that generate Co(i),n = ^75 for all n. Then, 
we obtain the fidelity of the QSE as 


F — + X] Crl + lC ")’ 
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For | a 2 1 1, it can be approximated as 
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